function [b , se , p , Radj ,T] = olsnan_gmm_EWC(LHS , RHS , c)
    T = nonnanrow_num([LHS , RHS]);
    nwlags = round(1.3 * T^(1/2));
    
    [b , ~ , ~ , Radj] = olsgmmnan(LHS , RHS , c , nwlags , 1);
    [~ , se , p , ~] = olsnan_EWC(LHS , RHS , c);%Calculate s.d. using EWC
end